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,^ ' The phases of primordial gravity waves is analysed in detail within a quantum me- 

chanical context following the formalism developed by Grishchuk and Sidorov. It 

<D '. 

is found that for physically relevant wavelengths both the phase of each individual 

mode and the phase difference between modes are randomly distributed. The phase 

X" 

sum between modes with oppositely directed wave-vectors, however, is not random 
and takes on a definite value with no rms fluctuation. The conventional point of view 
that primordial gravity waves appear after inflation as a classical, random stochastic 
background is also addressed. 



It is well known that gravitational waves (GW) can be produced in inflationary cosmolo- 
gies |IJ. They arise from quantum fluctuations in the gravitaional field on a curved de Sitter 
space-time. The generated perturbations are then often treated as a Gaussian noise with 
randomly distributed phases. During de Sitter expansion, it is believed that these fluctua- 
tions are then stretched "out" of the horizon where they "freeze" and remain constant until 
much later when they re-enter the horizon and appear as classical GW. These GW then 
form a classical stochastic background which fills the universe. It has, however, also been 
argued that the quantum-mechanical generation of perturbations, in constrast to a classi- 
cal, stochastic perturbation spectrum, actually possesses very specific statistical properties. 
Specifically, it has been argued that the phases of all modes of perturbations are essentially 
constant and fixed. Moreover, strong quantum correlations cause classical standing GW to 
form leading to important cosmological consequences ||. 

In this paper we shall present a detailed analysis of the statistical properties of the 
phases of relic GW from a quantum mechanical viewpoint. As the generation of these GW 
are fundamentally quantum mechanical in nature we shall thus be able to comment on 
the view that these GW subsequently form a classical, stochastic background. Our results 
show that the phase of each individual mode is random. More importantly, in all but one 
case, there is also the absence of any correlation between the phases of two different modes. 
Namely, both the phase sum and phase difference between any two aribitrary modes vanish. 
There exists, however, one exceptional case, that of two modes with oppositely directed 
wave-vectors, for which the phase sum of the two modes are highly-correlated. This is 
a manifestation of phase locking between modes which also occurs in quantum optics || . 
The phase difference between these two modes are still completely uncorrelated and random, 
however. In many instances, therefore, the relic GW may be approximated as a "classical" 
stochastic background. Nevertheless, the phase-sum locking asserts that the background is 
fundamentally quantum mechanical in nature. Based on these results, we shall comment on 
the usual classical treatment of relic GW. 

In our analysis we shall follow the framework developed by Grishchuk and Sidorov [|J. 
But instead of using 'standing wave' operators which decompose the two- mode squeeze 
operator into the product of two single mode operators, we shall stay with the original 



'traveling wave' operators where the physics is more straightforwardly seen. We warn the 
reader, however, that while the number operator for each mode is invariant under either 
choice, the phase operator is not and our definition of the phase operator differs from that 
given in |4j] 

The quantized graviton field operator hij may be written as 

ha^i) =C£IX® \^(v)e 1 ^ + ai\ V )e-^} , (1) 
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where C is an overall constant whose value is not imporant for our purposes, pfj(q) are the 
polarization tensors and s labels the two polarization states of the wave. a s drj) and aV(j]) 
are raising and lowering operators in the Heisenberg representation and with their evolution 
governed by the hamiltonian 

H = Y,i: {K-^K-W + Z*U^) a '-M + 2a ^) [a?fa)aVfa) - 4(77)0^(77)] } (2) 

q s=l 

where 17(77) = iBl /2R. Here R is the cosmic scale factor and the prime denotes derivative 
with respect to 77. 

Notice that the graviton field is coupled to the background metric, and thus to the ex- 
pansion of the universe, through a conformal time varying quadratic interacting hamiltonian. 
This hamitonian may be diagonalized and the Heisenberg evolution equations 

idaydrj = [at, H] , idap/drj = -[a^, H] , (3) 

solved exactly using a Bogolubov transformation. In fact, Grishchuk and Sidorov make the 
time dependent Bogolubov transformation: a%(r)) — > u q (r)) a~(r) ) + v q (rj)a s _^(rjQ), a s ^(r]) — > 
Uq(v) a P (Vo) + ^g( ? 7) a -g( ? 7o) where 770 is some initial conformal time. Requiring this to be a 
canonical transformation restricts 

Wt7)| 2 -K(77)| 2 = 1, (4) 

for all 77. Then, noting that in the presence of inflation R'/R — > as 77 — > —00, they choose 
as the initial conditions for the Bogolubov coefficients to be u q (r] ) = 1 and v q (r)o) = where 
7/0 is some initial time in the far past. 

While this is certainly a valid approximation, one eventually finds that it is more useful 

to write 



u q {rj) = e ie " coshr q , v q (r]) = e i{ - tq 2ct>q) sinhr g , (5) 

which explicitly satisfies eq. (f|). The functions r q {rj), e q (r]) and <j) q {rj) are called the squeeze 
parameter, rotation angle, and squeeze angle, respectively. If one then uses Grishchuk's 
initial conditions, one finds that r q (r)o) = and £5(770) = 0. The initial condition for <j> q is 
undetermineable, however, and it is this phase which will play a crucial role in determining 
the phase of the GW. A more careful analysis must there be done to establish the correct 
initial condition for it. 

To do so, let us begin by defining 

<x%n) = u q (v)a s $(Va)+v q (r])a-?(Vo) , 
ap(r]) = u q (r])ap (tj ) + 7^(77)^(770) , (6) 

where a|(?7o) and ay(r]o) are evaluated at some initial time 770. Then from the Heisenberg 
evolution equations, u q and v q evolve as 



du q 

dt] 


R' 

-iqilq + —V q , 


dVq 

drj 


R' 

~ lQVq + ~R Uq ' 



(7) 

The initial conditions for u q and v q are given at r/o and are fixed by requiring a^(rjo) and 
acp(r)o) to diagonalize H at this time: 

H(vo) = E E ^(<l) k%o)c4M + oi/(i»)ai^)l , (8) 
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since it is only in this way that a ground state for the system can be defined. From the 
standard Bogolubov analysis, we find that 



u q (Vo) = ^e W J 1 -^ , Vq { m ) = -e»J^, (9) 



and 

e{q) 



x 1 - f^H) ■ HO) 



Here 9 is an arbitrary phase for each mode q. We fix it by requiring that u q (r]o) — > 1 if the 
interaction term is turned off: |cr| — ^ so that 9 = — %/2 for all q. Notice also that because 



e(q) must be real number, the initial time 770 must be chosen such that 2 1 cr ( 770 ) | < q for all q. 
Usually tjq is chosen to be at such an early time that |cr(r/o)| ~ and this is not a problem. 



From Eq. (9), cosh 7^(770) = J(l + e)/2e, and e g (?7o) = 0. These argee with Grishchuk's 
initial conditions in the limit e(q) — ■> 1. The initial condition for 4> q (r]o) = — 7r/4 can now be 
determined, however. 

As was first noticed by ||, the above system is equivalent to what are called squeezed 
states in quantum optics. Namely, one notices that the tranformed operators are related to 
the original operators through a unitary transformation: 

ap(v) = ni( V )Si( V )ap( Vo )Si\ V )nf( V ) , (11) 

where 

Kftj) = exp {-ie q ( V ) (ap( Vo )ai( Vo ) + al/M^o))} , (12) 

is the two mode rotation operator while 

Sfa) = exp {r q ( V ) ( e -^Wa^o)a 8 » - e 2i ^ap( Vo )ai J( Vo )) } (13) 

is the two mode squeeze operator. 

It is currently held that quantum fluctuations in the gravitational field will be amplified 
as the universe undergoes a fast expansion. (See, for example || for a complete description 
of particle creation in the universe). This can be seen explicitly by looking at 

(0|r4|0> = (0|of fa)ajfa)|0) = \v q ( V )\ 2 = sinh 2 r g ( V ) , (14) 

for any mode q. It is then argued that at the present time these primordial, quantum 
fluctuations will appear as classical GW with a classical stochastic destribution of random 
phases. This, however, has never been explicitly established. 

It is known that the number operator and the phase operator 

exp{^|(ri)} = (/ + ntfrij) ' a|(ri) , exp{-^(r/)} = a| f ( V ) (/ + n|(^)) ' , (15) 

for the graviton field do not commute. The phase itself is not an hermitian operator, 
however, and suffers from a problem with multiplicity. We shall instead have to work with 
the operators 



cos<fl(r]) = - (exp{i^|(r?)} + exp{-z^|(^)} 

sinpffa) = — (exp{i(p a n (ri)} - exp{-z>|(7/)}) , (16) 

which are hermitian, and thus physical observables. One then finds that 

[ra|, cos tpQ = -i sin p| , [nf , sin y?|] = i cos y?| , (17) 

with the corresponding uncertainty relations 

An|ecosy>§ > -|(sinp|)| , Anfesin^f > -|(cosp|)| , (18) 



where A A = J (A 2 ) — (A) 2 is the rms fluctuation of the operator. The classical limit for 
massless particles, therefore, involves not only looking at the intensity of the radiation, but 
also at its phase. 

Using eq. (||), the rms fluctuation in the number operator for each mode is 

An| = \u k \\v k \ = V2coshr k (r))smhr k (i]) (19) 

so that An|/(n|) = v^2/tanhr fc . Now, it is usually stated |4| that although (0|/iy|0) = 0, 
because (0|n||0) 7^ 0, for k mode with a large squeeze parameter one can still view each 
mode of the GW as a classical wave with amplitude A^ = (OlnUO) 1 ^ 2 ■ One can certainly 
make this interpretation, but it does give the impression that the classical wave will have 
a well defined, definite amplitude. From eq. ([T9|), we see that this need not be the case. 
In fact, in inflationary cosmologies ||, modes just entering the horizon at the present day 
would have a r k ~ 120 giving Ant ~ (n|). The fluctuation of the amplitude of this classical 
wave will be on the order of its amplitude itself. In this case, we cannot characterize the 
resultant 'classical' wave as having any definite amplitude. 

Let us now consider the phases of these GW. We begin by calculating the average phase 
of each mode 

<0| COS ¥7jfo)|0) = (017^)^(77) 008^0)^(^)^(^)10) , (20) 

where we have used eq. (0). Since TZp\0) = |0), and by using the following factorization 

0, 



S} = — ^- exp {-a s J( Vo )a s _J( Vo )e^ tanhr,} 

exp {- a s ^(rio)a s ? (r]o) + ai/(»7b) a -$07o)] In cosh r 9 } 

exp {a^oK^oK 2 ^ tanhr,} . (21) 



Then 



<S| + |0> = — — £(e 2 ^ tanhr g )>,n> , (22) 

cosn r q n=0 

so that (0| cos ipUfl) |0) = 0. Similarly, (0| sin<^|(ry)|0) = 0. The average phase of any one 
mode vanishes, as expected for randomly distributed phases. This, however, does not tell 
us whether or not this random distribution is "classical". To do so, we calculate the rms 
fluctuation in the phase 

(Acos^) 2 = (0|(cos^) 2 |0) = \ - 1— |j- , (23) 

with (0|(cosy4) 2 |0) = (0|(sin^|) 2 |0). 

In eq. fl23|) we can see a deviation from the classical stochastic behavior. Suppose that 
the phase of this mode is classical in nature and can be described by a random, stochastic 
behavior. Denoting the phase of this mode by 9 q , then because everything are c- numbers, 
cos 2 9 q + sin 2 9 q = 1. If the phase of this mode is random, we would expect the stochastic av- 
erage (cos 2 9 q ) s to = (sin 2 9 q ) sto - Consequently, for a classical, random stochastic distribution 
of phase, one would expect (cos 2 9 q ) sto = 1/2. 

Since, however, (cos<^|) 2 and (sin^?.) 2 are operators, their sum need not add up to unity. 
As such, for a random distribution their expectation value (average) need not be 1/2, as it 
is for the classical stochastic distribution. In fact, any deviation from 1/2 is a sign of the 
quantum nature of the mode. As we can see from eq. (|23|), the quantum nature of the mode 
is always present, but becomes progressively smaller for large r p . Modes just entering the 
horizon at the present day have ar,~ 120 and for these modes (0|(cos<£>|) 2 |0) ~ 1/2. They 
are therefore essentially classical in nature and using a random, stochastic distribution to 
describe their phase would be correct. Notice, however, that for r q ~ 1, deviation from the 
classical behavior becomes pronounced and these modes are essentially quantum mechanical 
in nature. 



Of course, the absolute phase of any one mode is irrelevant. What is more interest is the 
relative phases between modes. Given the phase operator for any one mode, we follow |3j 
and define the phase sum and difference operators between any two modes as 

sin(<f s p ± (ft) = sin (ft cos f\ ± cos <^sin (f\, 

cos(y?|> ± (f\) = cos (p s f cos (pi if sin (ppsin (p\. (24) 

Notice that in the limit p — > q, t = s the sine difference operator does not vanish since 
[cos ip*, sin (pU 7^ 0. This once again underscores the fact that we are dealing with operators 
and not functions. 

It is then straightforward to show that for all t, s, and p ^ —q, 

(0| sin(^± <4)|0) = , <0| cos(<^± <p\)\0) = , (25) 

as expected since the average phase of each mode vanishes. What is of more interest is the 
expectation value of the squares of these operators 

(0|(sin(^-^)) 2 |0) = (0|(cos(^-^)) 2 |0> = i(tanh 2 r p + tanh 2 



(0|(sin(^+4)) 2 |0) = (0|(cos(^+4)) 2 |0) = ^(l + tanh 2 r p tanh 2 r, 
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(26) 



Once again in the limit of large r p and r q , these results go over to what one expects for a 
classical stochastic distribution of the phases. Consequently, we find that as long as p ^ —q 
the phase difference and sum between modes are completely random and in the limit of large 
r p and r q , can be accurately approximated as classical stochastic distribution of the phases. 
When t = s and p = —q the situation changes quite dramatically, however. First, we 
find that 

(0| sim>;— <^)|0) = 0, <0| cos^-v^lO) = 0, (27) 

as expected. But now 

(0| sm((fp + (ftp) |0) = sin 2cj) p tanhr p , 
(0| cos(y^ + (f s _p\ 0) = cos 2(j) p tanh r p . (28) 



Moreover, 



(sin(^— </- P -)r|0) = (0|(cos(^-^)) 2 |0> = -tanh 2 r p , (29) 

while 



{0\(sm(<pp+ <p a _ p )) 2 \Q) = - (l-tanh 2 r p + 4tanh 2 r p sin 2 (20 ? 



(0|(cos(<^+<^ p -.) 2 |0) = - (l -tanh 2 r p + 4tanh 2 r p cos 2 (20 p )) . (30) 

Once again, in the limit of large r p , we find that (Asin((pt— (p s _p)) 2 = 1/2, meaning that 
phase difference between the two modes is completely random. The fluctuation in the phase 
sum A sin(yj|.+ <£>i~-) = 0, however, and the two modes have a definite phase sum. This is a 
manifestation of two mode phase locking first calculated by [[| using distribution function 
methods. 

We thus see that although the phase differences between modes are uncorrelated and 
completely random, the phase sum between the p and —p modes with the same polarization 
is highly correlated with an average value of 2<p p and essentially no fluctuation whatsoever 
at large r p . This correlation was also founded by Grishchuk and Sidorov and was interpreted 
by them as the formation of a classical standing wave. Note, however, that this phenomenon 
is inherently quantum mechanical in nature and it is impossible to explain both eqs. ( |29| ) 
and (|30|) using classical stochastic arguments for the following reasons. 

Suppose that we wish to explain the results of the phase sum and difference analysis 
for p = —q in the large r p limit using classical stochastic arguments. Then eq. (|28|) implies 
that the phase distribution of both the p and — p modes are peaked at <ft p . That there is 
a width to this distribution can be seen in eq. (p9|) and in fact we see that there must be 
a randomly distributed background noise below <ft p . If, however, this background noise is 
present, then one would not expect eq. (|30|) to hold, since although we would expect the 
average phase sum to be 20 p , we still expect the noise to be present and would not expect 
the rms fluctuation in the phase sum to vanish. Classically, it should be 1/2 once more. 
Consequently, the results of the phase sum and difference analysis cannot be explained using 
classical methods. 

The phase locking of the p and —p modes also suggests the following method to measure 
the thermal history of the universe. Notice that for large r p , which correspond to most 
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physically relevant modes, the phase sum of the two modes is 20 p with essential no fluctuation 
whatsoever. From eq. (|7|), we see that the value of (f> p depends explicitly on the scale factor 
R which is in turn determined by the thermal history of the universe. Consequently, we 
propose to fix a specific direction along the celestial sphere and measure the phase sum of 
primordial GW with various momentum. This will determine each <p p . From eq. (0), r p can 
be infered by measuring the amplitude of the wave, although from eq. (|19D the fluctuation in 
this amplitude is expected to be large. Through these two measurements the thermal history 
along that direction can be extracted. Next, fix the magnitude of p and measure the phase 
sum = 20 p for a series of angles on the celestial sphere. If the universe was truly isotropic 
throughout its history, then p should be independent of angle. If not, then variation in (p p 
will be a measure of the anisotropy of the universe. In either case, a complete determination 
of <j) p and r p will be useful for probing the very early universe. 

How these experiments would be done is not as yet known. Not only hasn't any GW been 
detected yet (not to mention a primordial one), this experiment would also require measuring 
the quantum mechanical phase sum of two modes. It has only been very recently that the 
phase operator for optical waves has been measured experimentally ||, and no experimental 
measurement of phase locking has yet been found even for optical waves. Consequently, we 
would expect such experiments, if they can be done at all, can only be accomplished in 
the far, far future. Nevertheless, they have the potential to provide a very clean and direct 
measurement of the history of the universe. 

To conclude, we have in this paper completed a detailed quantum mechanical analysis of 
the properties of GW arising from fluctuations in the de Sitter vaccuum of the universe. We 
find that the present day charactorization of these waves as classical GW with a classical, 
stochastic distribution in the phase to be somewhat naive. First, while one can identify 
the rms fluctuation of each mode of the graviton field h%j as a classical amplitude, the 
fluctuation of this 'amplitude' is very large for large r p , the most physically relevant ones. 
Next, although the average phase of each mode seperately is randomly distributed and can 
be well approximated as a classical stochastic distribution, there are non-classical, strong 
correlations between the phases of different p. There is a definite value of the phase sum of 
the p and —p modes, although the phase difference of these modes, as it is between any two 
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modes, is completely random and stochastically distributed. This has been interpreted by 
Grishchuk and Sidorov as being due to the formation of standing waves in the universe. 

We now can understand the limitations of the usual classical treatment of the relic GW. 
In this treatment, the time evolution of the amplitude of classical GW generated from 
inflation is obtained from a two-point statistical average of a stochastic ensemble of classical 
gravitational fields which is then matched to the quantum mechanical two-point function of 
the gravitational field 0. The subsequent evolution of the classical waves is then described 
by the classical wave propagation in the expanding universe. This constitutes the so-called 
stochastic classical background of relic GW with randomly distributed phases. Intuitively, 
we say that during inflationary expansion fluctuations in the gravitational field will be red- 
shifted out of the horizon, after which they freeze and remain at a constant amplitude. Much 
later when they re-enter the horizon during either the radiation- or matter-dominated era, 
they will appear as classical oscillating GW. 

There are, however, two different phases which must be considered if one wishes to view 
the system in this way. The first is the temporal phase of oscillation of waves. This is 
precisely fixed in inflationary cosmologies. Once a frozen mode starts to oscillating as a 
classical wave at the time it crosses the horizon and re-enters the universe, the temporal 
phase is completely determined by the classical equation of motion. However, inflation does 
not predict the location of the nodes of these frozen modes since quantum fluctuations assign 
equal probability to modes which differ only by a spatial translation. Herein lies the other 
phase which must be considered: that due to the spatial part of the GW. It is this phase 
which is completely random. Consequently, in this classical treatment the overall phase of 
GW are randomly distributed. 

Based on our results of the statistical properties of the phases of GW, the above classical 
description has to be amended. Since the phases of different modes are almost uncorrelated, 
in many aspects the relic GW can be well as a stochastic classical background radiation. 
Phase-sum locking, however, is still present and is a strong indication of its inherent quantum 
character of these GW. Moreover, this seperation of the phase of each mode into spatial and 
temporal parts is completely artificial and araises only from the desire to use a classical 
description and interpretation of the quantum two-point function. From the experience of 
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calculating the vacuum expectation values of the phase operators, we can hardly separate 
the temporal phase from the spatial phase. This is not surprising, though, since there is 
only one single overall phase in the quantum mechanical treatment. Consequently, although 
the classical approach advocated in is quite physical and is, as we have seen, a good 
approximation of an inherently quantum mechanical phenomenon in many instances, it can 
not fully describe the properties of the relic GW. Our results suggest that one should instead 
use the quantum mechanical approach to deal with these GW. 
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